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ABSTRACT 


In this report, a phase measurement is added to the Underwater Sound Reference Division 
(USRD) reciprocity coupler primary calibration procedure for an H48 reference 
hydrophone. Data acquisition equipment is added to record time-series data from the 
hydrophone under test and from the reciprocal transducers. The complex-valued 
hydrophone sensitivity is calculated. The sensitivity magnitude is compared to 
measurements from the standard coupler calibration procedure, and the complex sensitivity 
data are also fitted to a simple high-pass circuit model. The model is used to estimate the 
low-frequency cutoff of H48 hydrophone SN4. The low-frequency cutoff measured in this 
report is about 0.2 Hz higher than that originally measured and specified when the H48 
hydrophones were first built. The new results show significant roll-off in phase below 
10-20 Hz, a range where the phase is typically assumed flat during the standard calibration. 
By 1 Hz the phase roll-off is about 20°. The error analysis of the original coupler is 
summarized and error and uncertainty due to new data acquisition equipment and phase 
measurement added. Some errors due to simplifications in the acoustics of the coupler are 


left to future work. 
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CHAPTER 1: 
Introduction 





The Underwater Sound Reference Division (USRD) reciprocity coupler is used to cali- 
brate the reference hydrophones that disseminate the primary calibration to the secondary 


references used by the U.S. Navy (USN) and other organizations. 


To completely reproduce an acoustic waveform during calibration, the magnitude and phase 
of the waveform are required. Historically, hydrophone calibration is more often performed 
gathering only magnitude information from the hydrophone under test. Phase response is 
more difficult to measure accurately and can be assumed constant when hydrophones are 
operating well below their resonance frequencies and above their low-frequency cutoff. The 
current USRD coupler calibration procedure for the H48 reference hydrophone measures 
magnitude only and assumes a constant phase response over the bandwidth from 5 Hz-2000 
Hz. The reference hydrophones are operating well below their first resonance, but the bottom 
of the frequency range can be affected by the hydrophone’s low-frequency cutoff, especially 
for tests where secondary references are used to measure frequencies below 5 Hz. Adding 
a phase measurement to the calibration would reduce the phase error at the low end of the 
calibration frequency range where the phase response is no longer constant. These phase 
errors can undermine accurate direction finding and degrade broadband measurements made 


by sonar arrays. 


As demand rises for phase measurements, an accurate phase calibration of USRD primary 
references becomes more important. The USN can use this information to enhance ac- 
curacy in beamforming, in broadband applications, and in low-frequency sound pressure 
measurements required by the Office of Naval Intelligence (ONT). Modern towed arrays 
could be supplied with more accurate calibrations at the low end of their operating frequency 
band. Phase information will allow for measurement of the resistor-capacitor (RC) time 
constant of the equivalent single-pole high-pass filter circuit at the input to the preamplifier 
in a hydrophone channel of a towed array, an important quality control parameter. An 
out-of-specification RC time constant can lead to phase response errors between otherwise 


identical hydrophones. This constant increases as parasitic resistances in the hydrophone 


increase and as the input resistance to the preamplifier decreases. An RC time constant 
that increases over the time can render an array unable to accurately measure signals it was 
previously calibrated to measure. Phase calibration of hydrophones will complement the 
magnitude calibration and allow the RC time constant to be measured and tracked without 


taking the array apart. 


Hydrophone phase calibrations can be performed in a variety of ways. Luker and Van 
Buren [1] performed a phase calibration in an open tank using a three-transducer reciprocity 
method. The projector, hydrophone, and reciprocal transducer were all arranged in a straight 
line, with the hydrophone in the center. The hydrophone could be rotated in place to face 
either the projector or reciprocal transducer or be removed from the tank, eliminating phase 
errors due to distance between the transducers and hydrophone, and sound speed. Hayman 
et al. [2] performed a free-field reciprocity calibration, with phase, in an open tank facility, 
then compared the results to a phase calibration using optical interferometry. In that setup, 
a thin, acoustically transparent, optically reflective membrane called a pellicle is placed a 
fixed distance from a source and an optical interferometer or laser vibrometer is used to 
measure acoustic particle velocity at the pellicle. The pellicle is then replaced with the 
hydrophone under test, and hydrophone voltage measurements taken with the same source 
are used with the acoustic pressure to calculate sensitivity. The location of the acoustic 
center of the hydrophone must be carefully matched to the location of the pellicle. 


The USRD reciprocity coupler provides an accurate and consistent primary calibration up to 
very high pressure that is not subject to many of the sources of error in the above mentioned 
methods, due to the equipment setup and the physics of the coupler procedure. This paper 
explains how a phase measurement was added to the USRD reciprocity coupler. In short, 
the magnitude and phase of all electrical parameters in the reciprocity calibration procedure 


must be measured and the complex-valued reciprocity parameter must be used. 





CHAP TER: 
Theory 





The reciprocity coupler is a fluid-filled cavity in which three transducers are installed to 
perform a reciprocity calibration to calculate the low-frequency, free-field sensitivity of an 
unknown hydrophone. Figure 2.1 shows the cross section of an example coupler. Identical 
transducers are installed on the left and right sides of the coupler and a reference hydrophone 
(to be calibrated) in the center. Castor oil fills the sealed cavity. 
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Figure 2.1. A cross-section view of three transducers installed in an example 
reciprocity coupler. Source: [3]. 
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The pressure sensitivity is obtained via the reciprocity method by measuring the current to 
the projector and the voltages at the receiving transducer and unknown hydrophone. The 
reciprocity parameter ties these measurements to the sensitivity of the unknown hydrophone. 


The coupler is designed around a set of assumptions that, when satisfied, ensure the pres- 


sure and free-field sensitivities are equal, allowing the more useful free-field sensitivity 


measurement to be calculated during the calibration. 


Throughout this report, complex quantities are shown in bold font. 


2.1 Coupler Calibration Assumptions 
Several assumptions are made when performing a calibration in a reciprocity coupler to 
enforce the condition that the sound pressure is the same everywhere in the coupler at a 


given instant in time, as described in [4]: 


1. The wavelength of sound must be much larger than the largest dimension of coupler, 
its length. 

2. The walls of the coupler and transducer elements must be rigid — they must have a 
much higher impedance than the acoustic medium. 


3. No air or any other low impedance material can be present in the coupler. 


The hydrophone sensitivity measured in the coupler is the pressure sensitivity, but the free- 
field sensitivity is the desired measurement. Fortunately, in the conditions of the coupler, the 
pressure sensitivity and free-field sensitivity are equivalent. Figure 2.2 shows a hydrophone 
with an open circuit voltage, e, and an acoustical impedance, Z,. The applied pressure 
at the acoustic input in the coupler is pg and an acoustical Thevenin generator represents 
the effects of a free field: the radiation impedance, Z;, and the blocked and free-field 
pressures, pp and pr. The ratio a is called the diffraction constant, D, although it varies 


with frequency. 





Acoustical Thevenin Generator 


Figure 2.2. An equivalent circuit of a hydrophone in a free field. Source: [4]. 


In order for p ¢ and py» to be equal, the diffraction constant D must be equal to 1, a condition 
satisfied if the hydrophone is much smaller than an acoustic wavelength, which must be true 
if the wavelength of the sound is much larger than the coupler cavity. Analyzing the circuit 
diagram, py» and pq can be considered equal if Z; is much smaller than Z,. In situations 
where the hydrophone’s physical dimensions are much smaller than the wavelengths, Z, 
will be small and the Z, << Z, criteria will be met by a piezoelectric transducer operated 
at frequencies away from resonance. With these two criteria met, the free-field pressure is 


equivalent to the applied pressure at the acoustic input [4]. 


2.2 Reciprocity Calculation for the Complex Voltage Sen- 
sitivity 
The primary calibration performed using the USRD reciprocity coupler takes advantage 
of the reciprocity method to calculate the complex voltage sensitivity of an unknown 
hydrophone. In a reciprocity calibration, the voltage sensitivity of a hydrophone can be 
calculated through the measurement of six electrical parameters made on three transducers 
(reciprocal transducers X and Y, anda hydrophone) as shown in Table 2.1. Two measurement 
setups are required to obtain all six electrical parameters, (a) reciprocal transducer X 
transmitting and reciprocal transducer Y and hydrophone receiving, and (b) reciprocal 
transducer Y transmitting and reciprocal transducer X and hydrophone receiving. The 
electrical parameters are: ix, the input current to transducer X when it is transmitting, eyy, 
the open-circuit output voltage from the hydrophone when transducer X is transmitting, 
eyx, the open-circuit output voltage from transducer Y when transducer X is transmitting, 
iy, the input current to transducer Y when it is transmitting, eyy, the open-circuit output 
voltage from the hydrophone when transducer Y is transmitting, and eyy, the open-circuit 


output voltage from transducer X when transducer Y is transmitting. 


The output voltage measurements can be expressed in terms of a voltage sensitivity and the 
pressure transmitted by the transmitting transducer. They are manipulated with the goal 
of expressing the sensitivity of the unknown hydrophone while eliminating any acoustical 
parameters of the other transducers in order to produce a primary calibration as shown 


in [4]. For configuration (a) this expression is 


enx = Mupx = MnixSyx, (2.1) 
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Table 2.1. Reciprocity calibration measurement groupings. Adapted from [4]. 








. Input Current Transmitting Receiving Output 
Configuration voltace 
Measurement Transducer Transducers 
Measurements 
oe Repl gmc’ enn 
2 Transducer X P HX, €yx 
Transducer Y 
(b) i Reciprocal es es : , 
: Transducer Y P HY, €XY 
Transducer X 





where My is the open-circuit receiving voltage sensitivity of the hydrophone and Sy is the 
transmitting current response of reciprocal transducer X. Unlike the free-field form of this 
expression, there is no dependence on separation distances, and the transmitting current 
response is the pressure inside the cavity for a given current — there is no reference distance 


term. Similarly, the open-circuit output voltage of the transducer Y is 
eyx = MyixSx, (2.2) 


where My is the open-circuit receiving voltage sensitivity of reciprocal transducer Y. The 


ratio of the output voltages in configuration (a) is 


eHx Mu 
— 2.3) 
eyx My : 
From configuration (b) 
eny = MniySy , (2.4) 


where Sy is the transmitting current response of reciprocal transducer Y. Using the property 
of the reciprocal transducer that the ratio of the receiving sensitivity and transmitting 


response is equal to the reciprocity parameter, 


= 3 (2.5) 


J 


which is only a function of acoustical properties of the coupler fluid medium. Equations 


2.3, 2.4, and 2.5 can be solved for the hydrophone sensitivity as shown in Equation 2.6. 


e e 
My = ,|——** 7. (2.6) 
eyxly 


In the USRD coupler calibration procedure, redundant measurements using reciprocal 
transducers X and Y as both a projector and receiver reduce measurement uncertainty and 
create convenient symmetry to simplify phase measurement. Measurement (b) in Table 
2.1 is collected when the cables driving the projecting transducer and recording from 
the receiving transducer are interchanged from configuration (a), and the measurement is 
repeated. The extra output voltage measurement, e yy, is incorporated into the denominator 


of Equation 2.6 where the geometric mean of both receiving transducer measurements gives 


e e 
My = ,|—2*—* sy, (2.7) 
Veyxlx X exyly 


In this coupler setup, current measurements are obtained by measuring the voltage across 
a standard capacitor placed in series with the input current. Using the relationships iy = 
€cxjwCp, and iy = ecyjwC,, where ecy and ecy are the voltages across the standard 
capacitor and C, is the capacitance of the standard capacitor, the currents in Equation 2.7 


can be substituted for voltages, producing 






CHXeHY J 


Veyxecx X exyecy JwCy 





(2.8) 


2.3 Reciprocity Parameter 

The complex pressure reciprocity parameter, J, must be calculated for a small rigid- 
walled cavity in order to obtain the hydrophone sensitivity from the voltage and current 
measurements in Equation 2.8. The beauty of the reciprocity parameter is that it depends 
only on acoustic properties of the fluid medium. No electrical properties of any of the 
transducers have to be known, making the reciprocity method a primary calibration. To 
derive this parameter, the definition of the reciprocity parameter as the reciprocal of the 


acoustical transfer impedance [5] is used. 


First, assume a simplified coupler geometry as shown in Figure 2.3. This cross section of 
a right cylindrical cavity of length L is capped by identical, reciprocal transducers at either 
end. Assume plane wave motion along the x axis with no motion in the radial or angular 


directions. The acoustic field can be modeled by the network diagram in Figure 2.4. 





Figure 2.3. A simplified reciprocity coupler geometry consisting of a cylin- 
drical cavity capped by two reciprocal transducers shown in cross section. 





Figure 2.4. A network diagram of acoustic field inside a cylindrical cavity 
capped with two reciprocal transducers. 


The pressure and volume velocity can be related by the mutual acoustical radiation 
impedance matrix in the network equations 


Pi U; 
= : 2.9 
e (i | as 


In this matrix, Z;; is the ratio of the acoustic pressure at transducer i to the volume velocity 


Zi, Zi2 
Zo, Zr 








of transducer j, with all other transducers but transducer j blocked. Z1; and Z22 are "self" 
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radiation impedances, the acoustical radiation impedance seen by one transducer when it 
alone is radiating; and Z)2 and Z2 are "mutual" or "transfer" radiation impedances, the 
ratio of the blocked pressure at one transducer to the volume velocity of the other. 


To derive expressions for the terms in the acoustical radiation impedance matrix, start with 


the relationship for the pressure of a plane wave, 
p(x,t) = (Acos(kx) + B sin(kx))e/, (2.10) 
where A and B are constants, k is the wave number, w is angular frequency, and f is time. 


From there, the particle velocity is given by 


1 0 k : 
u(x,t) = pe ae = ——(Asin(kx) — Bcos(kx))el. (2.11) 
Jwp ox jwp 


The network variables p;, po, U;, and U2 (Figure 2.4) are related to p(x, t) and u(x,t) by 


Pi = p(O,t), po = p(L,t), 
U, = Su(0,t), U> = —Su(L,t), 


(2.12) 


where S is the cross sectional area of the cavity. Because the acoustical radiation impedance 
matrix is symmetric, only one self- and transfer impedance needs to be evaluated, and the 
solution can be applied to the other. Focusing on the problem where there is a transducer at 


x = L and a rigid wall at x = 0, the boundary condition exists 
U@~=0)=0.. B=0. (2.13) 
Note that with B = 0, the relationship in Equation 2.10 describes a standing wave. 


Solving for the radiation impedance, 


Po ~=—~piL,t) = Acos(kL)el” jpc 
Z = Z SS SS SS SS SS ee t kL ; 2.14 
ee" Oy SUED “ESA sin(k L)e/" oor h) an 


and for the acoustical transfer impedance, 








P\ Pi P2 
Z = Z — — 
12 = 4 = 4 we 
= PLIPS cot(kL) 
p2 S 
0,1) j 
oj BAND TOE sce 
p(L,t) § 
Ae! J pec 
Se | cote 
Soe S } cot ) 
J pc 
= —- KL). 
5 csc(kL) 


Using the first two terms of a Taylor series expansion gives the following long-wavelength 


approximations for the self-radiation impedance 


Zo. = Z 4, = — 





Jpc (1 kL Baa jwpLl 
SHEN 2.15 
S (az 3 joSL” 3 ve) 


and the acoustical transfer impedance 





Zin = Zn = (zz + = wt ere (2.16) 
where Bug = pc? is the adiabatic bulk modulus of the fluid. The first term in each transfer 
impedance shown in Equations 2.15 and 2.16 contains the bulk modulus in the numerator, a 
term that represents the elasticity of the coupler volume under pressure, and the frequency 
and volume in the denominator. This term represents the acoustical compliance of the fluid. 
The second term in each impedance contains the density and length in the numerator and 
the surface area in the denominator, forming an inertance term that describes the pressure 


difference required to cause a change in flow rate with time. 


The reciprocity parameter, J, is equal to the reciprocal of the acoustical transfer impedance 
[5]. Using the first term of the Taylor series expansion in powers of kL for Z12, from 
Equation 2.16, the commonly accepted reciprocity parameter for the coupler is 


1 jwSL | 


J=- 
Zin pe? 


jwC, (2.17) 
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where C is the acoustical compliance of the coupler fluid medium as long as the fluid 
compliance is much greater than the compliance of the coupler walls and transducers. The 
coupler is compliance (or stiffness) dominated as long as a >> Mh For the worst case 
conditions at the highest temperature, lowest pressure, and highest frequency in the H48 


hydrophone test in this report, a = 1.32 and Mh = 0.126. 


From Equation 2.8, substituting the complex pressure reciprocity calibration parameter 


from Equation 2.17, we get 


V 
Mu = _ eee —_ ("| V/Pa, (2.18) 
Veyxecx X exyecy \Cppc 


where V is the volume of the coupler cavity. The voltage sensitivity in dB re 1Vju/Pa is 
given by 
My = 20log(|My|) — 120 dB re 1V/pPa. (2.19) 


The factor of 120 is used to convert the dB reference to uwPa. This voltage sensitivity is 


equivalent to the free-field voltage sensitivity while the assumptions in Section 2.1 are true. 


11 


THIS PAGE INTENTIONALLY LEFT BLANK 


12 





CHAP ERS: 
Experimental Setup and Methods 





A reciprocity calibration, including phase, was performed using the USRD H48 reciprocity 
coupler and supporting equipment. New data acquisition equipment were used to obtain 
the time series data necessary to add a phase measurement. A least-squares method was 
used to process the time series data to obtain magnitude and phase. Systematic errors 
in the experimental setup were calculated and removed and a simple, equivalent-circuit 
high-pass filter model was fitted to the data and used to estimate the cutoff frequency of the 
hydrophone under test. 


3.1 Physical Setup 

A diagram of the physical setup is shown in Figure 3.1, with component details provided in 
Table 3.1. A signal is generated at the function generator, amplified, and passed through a 
standard capacitor wired in series with the projector. A custom isolated differential amplifier 
measures the voltage across the standard capacitor, which is proportional to the current to 
the projector. A hand operated pressure ram controls the hydrostatic pressure in the coupler 
fluid, castor oil. Chamber temperature is controlled by heating or cooling distilled water 


from a temperature bath. 


Table 3.1. Equipment List 





Diagram # Component Make/Model 
1 HP 3325A 
Krohn Hite 7500 
Standard 1 yw F capacitor, custom differential amplifier 
USRD H48 Reciprocity Coupler 
Ithaco 1201 
NI BMC-2110 and NI-6259 
Laptop Computer with Matlab 
Pressure Ram 
Neslab Endocal RTE-8DD 





OMANI DANBRWN 
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3. Capacitor/ 
2. Amplifier Isolated 
Amplifier 


1. Function 
Generator 


5. Preamplifier 


6. Data 
Acquisition 





8. Pressure Analog Signal 
Ram 9. Temperature === Digital Signal 
Pea — fluid 


Figure 3.1. Block diagram of the experimental setup 


Within the coupler, three transducers take part in the H48 USRD calibration: two reciprocal 
lithium sulfate spheres and the H48 reference hydrophone designed for the coupler, which 
also has a spherical sensing element. The reciprocal transducers are mounted on threaded 
steel plugs and the housing of the primary reference is designed for the threaded hole in the 
top of the coupler, which positions the H48’s sensing element between the two reciprocal 
transducers. The receiving transducer is wired to a preamplifier after which the signal is 
recorded. The wiring to and from reciprocal transducers X and Y can be interchanged in 
order to record signals generated by either reciprocal transducer. The hydrophone has a built 
in preamplifier and is connected directly to the data acquisition equipment. A computer 


with Matlab is used to record the time series data and control the function generator. 


3.2 Complex Signal Measurement 
The magnitudes and phases of the voltages and currents required to calculate the complex 


sensitivity of the unknown reference hydrophone will be obtained using a least-squares 
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estimation process originally developed by Petr Vanicek [6]. Least-squares methods have 
been shown to perform better than conventional spectral analysis when the signal of interest 
is contaminated with systematic noise such as a direct current (DC) bias shift, a linear trend, 
or a sine wave of known frequency, even when that frequency is very close to that of the 


signal [7]. 


The measured voltage (or current) signals can be represented by 


e(tn) =C + Acos(2m ftn + d) 
= C + Acos(¢) cos(2aftn,) — Asin(¢) sin(27 ft,), (3.1) 


where e(t,) is the measured signal at any discrete time instant, t,, A is the amplitude of 
the signal, ¢ is the phase of the signal, and C is a DC offset. As shown in [8], in the 


least-squares estimation, the design matrix is built from the equation 


e(tn) = Bo + Bi cos(27 fty) + Bosin(27ftn), (3.2) 


where the beta coefficients calculated by the least-squares fit can be used to calculate the 
DC offset, signal amplitude, and signal phase by comparison with Equation 3.1. 


The beta coefficients can be used to reconstruct a complex signal, e, given by 


e = Ae!®, (3.3) 
where 
A= 4p? + 8 (3.4) 
and F 
-1 | ~P2 
= tan |——}]. 3:5 
( Bi ve 
The complex signal representation 
e= fi +jP2 (3.6) 


is used to calculate hydrophone sensitivity rather than magnitude and phase separately 


because it allows for more convenient processing and is necessary for proper accounting of 
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measurement uncertainty, discussed later. 


This model is used in the typical least squares form 


y = XB, (3.7) 
where 
e(to) 1 cos(2mfto) sin(2z fto) 
_ | ett) —s eos(2m fit) sin(2n ft) | da Be a | ae 
B2 


e(tn) 1 cos(2rftn) sin2afty) 


The vector of coefficients 6 can be solved for given a vector of values y with a matrix of 
values X as 
B= OCD AX 7. (3.9) 


3.3. Systematic Errors and Corrections 

The systematic errors in the coupler measurement apparatus contribute more to the total error 
than the uncertainties in the measurements and their relevant constants. These (systematic) 
errors arise from simplifying theoretical assumptions not fully borne out in the measurement 


setup and phase error introduced from the phase responses of the measurement equipment. 


Systematic phase corrections are applied to 2.18 as given by 


MPhCor — €nxeny ( V | ; 
8 SS a, 
V(Go-3ae yx) (Go-3necx) X (Gr-30e xy) (Go-3necy) \Cppc 


(3.10) 
where the phase shifts associated with complex gains G2_3, and G2_3, are defined in Section 


3.3.1. The magnitude of each of these complex gains is unity. 


Systematic magnitude errors are treated by Zalesak [9] and summarized in Sections 3.3.2 


and 3.3.3. Those corrections are applied after the phase corrections in Equation 3.10 
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producing 


MagCor _ Fa Foy lMiicer 
=.= H . 


V For 


F is applied to the compliance term V/pc? and is derived in Section 3.3.2; Foy is applied 


juror” (3.11) 








once to eyx and once to ey and is derived in 3.3.3; and Foy is applied once to eyx and 


once to exy and is also derived in 3.3.3. 


3.3.1 Phase corrections 

To isolate the phase response of the hydrophone, phase change contributions by the other 
components must be measured and eliminated. The block diagram in Figure 3.2 illustrates 
the coupler setup, which can be considered a linear electro-acoustic system. Significant 
measurement node locations are shown as red dots and the signal path as solid black 
lines. Dashed lines indicate the different measurement pathways through the system. Each 


component in the diagram has an associated linear complex gain element. 


Capacitor/ 
Amplifier Isolated 
Amplifier 


Function 
Generator 


Preamplifier 


Data 
Acquisition 
—— — —>(a) (b)<---- 





=—_——— Signal Path 
===: Measurement Path Direction 
e Measurement Nodes 


Figure 3.2. Reciprocity coupler system block diagram with measurement 
nodes indicated by red dots and phase changes indicated by dotted lines 
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Three measurement pathways through the electro-acoustic system are of interest: (a), (b), 
and (c). Each one begins at measurement node (1) in the fluid in the coupler cavity at the 
face of the hydrophone and ends at measurement node (4) where the hydrophone output 
voltage is recorded at the data acquisition system. Pathway (a) leads from the fluid in the 
coupler cavity through the receiving reciprocal transducer (Y in Figure 3.2), preamplifier, 
and data acquisition system. Pathway (b) leads from the fluid in the coupler cavity through 
the transmitting reciprocal transducer (X in Figure 3.2), capacitor/isolated amplifier, and 


data acquisition system. The complex gain through pathway (a) or (b) will satisfy 
Grot = Gi-2m * Go-3m + G3-4m, (3.12) 


where subscript number signify the nodes the complex gains act between and m can signify 


pathway (a) or (b). The gain magnitudes satisfy 


Gisim = Gi-tm* G2-3m * G3=4m (3.13) 


and phases satisfy 
Protm = $1-2m + $2-3m + $3-4m- (3.14) 


The change in phase in pathway (c), the unknown hydrophone phase response, is the phase 
of the complex sensitivity of the hydrophone calculated in Equation 2.18 where the signals 
recorded at measurement nodes (3a), (3b), and (4) are eyy, ecx, and ecy, respectively. 
The other parameters in Equation 2.18 are recorded separately when the measurements 
are repeated with the preamplifier connected to transducer X and the capacitor/isolated 
amplifier connected to transducer Y (not shown in Figure 3.2). The phase change, orm for 
pathways (a) and (b) must be calculated and applied as a systematic phase correction to the 
signal measured at nodes (3a) and (3b) before they are used in Equation 2.18. 


The following subsections discuss each of the complex gain terms. 
Gi-2m 
As the signal passes from measurement node (1) in the coupler fluid to measurement node 


(2a) or (2b), the voltage output and current input to transducer Y and X, respectively, it 


is modified by the complex gain represented by the receiving sensitivity and transmitting 
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current response of reciprocal transducers Y and X. These gain and phase changes will 
not affect the hydrophone sensitivity calculation because the reciprocity method allows 
the complex receiving sensitivities and transmitting current responses of the reciprocal 
transducers to be rewritten in terms of the reciprocity parameter, which is equal to purely 


acoustic properties of the coupler, as shown in Equations 2.1-2.6. 


G2-3m 

As the signal passes from measurement nodes (2a)/(2b), the electrical output/input to 
the reciprocal transducers, to (3a)/(3b), the outputs of signal amplification equipment, a 
complex gain is applied equal to the complex response of the amplifying equipment. The 
preamplifier and isolated amplifier are both unity gain, but they do shift the phase of the 
signals they amplify. To measure phase change ¢2_3,, the phase response of the preamplifier 
was measured. The phase response of the Ithaco preamplifier is shown in Figure 3.3 and is 
nearly flat within the +0.5°accuracy of the HP35565A used to make the measurement. 
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Figure 3.3. Ithaco 1201 frequency response 


The same procedure used to measure $2-3, is applied to the standard capacitor/differential 
amplifier package to measure phase change ¢2-3,. The standard capacitor is wired in 
series with the projector. The isolated amplifier measures the voltage across the current- 


measuring capacitor. The phase response of interest is that between the signal leaving the 
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capacitor for the projector and the amplified voltage signal measured across the capacitor 
at the isolated amplifier output. The package has three accessible ports: the input to the 
standard capacitor (C;,,), the output from the standard capacitor to the projector (C,,;), and 
the voltage measurement output from the isolated amplifier (V,,;). The response of interest 
(et) could not be measured directly with the available equipment, so two responses were 
measured: (1) injecting a signal into the capacitor input and measuring the signal at the 


amplifier voltage measurement output (Yau) and (2) injecting a signal into the capacitor 
Cou 
Cin 
the desired frequency response given by the equation below and shown in Figure 3.4. 


input and measuring the signal at the capacitor output ( 





). Dividing the two produces 











C Cout 

a ; 3.15 
Vout Vout ( ) 
Cin 
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Figure 3.4. Frequency response of standard capacitor/isolated amplifier 
package 
G3-4m 


As the signals pass from measurement nodes (3a)/(3b), each at a different channel of the 


data acquisition system, to node (4), they are subject to a complex gain that represents 
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magnitude and phase changes that occur due to unequal treatment of the signals at different 
channels. It is assumed that there is no change applied to the magnitude of the recorded 
signals from one channel to the next, but the same is not true for phase as the channels are not 
recorded simultaneously. Inspecting Equation 2.18, however, shows that any pair of equal 
and opposite phase adjustments (e? terms) made to signals from the reciprocal transducers 
in the denominator will cancel. The data acquisition system records samples from channels 
sequentially from the first to the last active channel, then repeats. With the hydrophone 
positioned in channel two, the voltage from the receiving reciprocal transducer in channel 
one, and the voltage from the current measuring capacitor in channel three, phase changes 
imparted by the data acquisition system will be made in equal and opposite amounts to the 
reciprocal transducer measurements in the denominator of Equation 2.18. Therefore 63-44 


and $34) cancel out. 


If the channels on the data acquisition box, an NI-6259, cannot be arranged symmetrically 
around the hydrophone input channel, the time difference between channel scans can be 
measured and the associated phase correction applied. The NI-6259 times acquisition using 
the sample clock for sample rate and the convert clock to time the sequential scans over the 
active channels within each sample. The driver software places 10 us of padding between 
channel scans unless there is not enough time between samples to pad between each active 
channel, at which point the channel scans are divided equally over the time period between 
samples. Using three active channels with the sample rate set to five times the signal 
frequency, the coupler measurements all receive 10 js padding between channels from 1 
Hz to 2 kHz. This padding has the desirable effect of eliminating measurement error due to 
settling time in the NI-6259, which (time) ranges from 1-2 us [10]. The padding time was 
measured by placing a T-fitting on the signal output of the function generator and recording 
the function generator signal in two adjacent channels of the NI-6259. Dividing the phase 
difference between the two recorded signals by the frequency gave an average time delay 
between channels of 11 us: 10 us padding and 1 us to record the measurement. Over 
the 5 Hz—2 kHz range of the H48 hydrophone, the associated delay between channels was 


constant. 
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3.3.2 Systematic error from compliant transducers and coupler walls 
The acoustic compliance found in Equation 2.17 was derived assuming that the transducers 
were rigid. In actuality, transducer compliance adds a small amount to the acoustic compli- 
ance term. The acoustic compliance of the transducers can be calculated using the change 
in outer radius for a thick walled hollow sphere for a change in external pressure, given 
by [11]: 
—qa (“ — v)(b* + 2a?) 
Aa = — |-——>.—__— — VI 


E 2(a3 — b3) ere) 


where Young’s Modulus, E, is 6.45 x 10!° N/m? for Navy Type I lead-zirconate-titanate, 
Poisson’s ratio, v, is 0.343, and the inner and outer radii, a and b, are 9.53 mm and 12.7 
mm, respectively. Using g = 1 Pa to calculate the change in radius per Pa, the compliance, 
or change in volume per Pa of the two reciprocal transducers and reference hydrophone are 
each 4.12x10!° m3/Pa, adding 1.236x 10-5 m3/Pa to the overall coupler compliance [9]. 


The compliance of the coupler wall, also assumed to be rigid, introduces error into the 
measurement. The coupler wall compliance was estimated using a finite element model [12] 
to be 2.66 x 107! m?/Pa with an estimated uncertainty of +2.7 x 107!°. A check on the 


finite element analysis can be calculated using Equations 3.17 and 3.18 from [11], 


_ qba(1t+yv)+b(1 —2y) 


A 
? E a2 — b* 


(3.17) 


and 
_ gl (1 - 2y)) 


E @-b ’ 


to determine the change in volume of a cylindrical chamber under internal pressure, assum- 


Al (3.18) 


ing the ends of the cylinder are rigid. 


The coupler chamber is made of type 304 stainless steel, with a Young’s modulus, EF, 
of 1.903 x 10!1 N/m? and Poisson’s ratio, v, of 0.305. The main cylindrical portion 
of the coupler is 86.77 mm long, /, with an inner radius, b, of 22.23 mm. Since the 
coupler has a rectangular cross section, a mean outer radius, a, of 61.32 mm is used to 
calculate an equivalent cylinder wall thickness. The compliance contribution of the main 
cylindrical portion, calculated as the difference in volume per unit of applied internal 


pressure (g = | Pa), is 2.252 x 10-!5 m3/Pa. The contribution of the cylindrical side arm 


pans 


that contains the hydrophone can be calculated by multiplying the main cylinder compliance 
by the ratio of the total volume to the main cylinder volume without transducers installed: 
1.530 x 10-4m? and 1.347 x 10~4m?, respectively. The total compliance, accounting for the 
hydrophone side arm is therefore 2.559 x 107!>m?/Pa, which agrees with the finite element 


analysis within the modeling uncertainty [9]. 


The compliance of the transducers and coupler walls add to the castor oil compliance to 
give the total compliance value that should be used to calculate the hydrophone sensitivity, 
5.643 x 10-'4m3/Pa. Equation 2.19 assumes only the castor oil has a compliance, so this 
addition is equivalent to multiplying Equation 2.8 by the factor 


2 

pc 
|e ean 3.19 
A ae (a1) 


where C7 is the compliance of the transducers and coupler walls. This correction amounts to 


between 0.25 and 0.33 dB over the temperature and pressure range of the H48 calibration [9]. 


3.3.3. Wavelength correction factor 

In developing the sensitivity equation for the coupler, it was assumed that pressure was 
constant everywhere in the coupler. This statement is a good approximation, but at higher 
frequencies there will be a non-negligible pressure variation across the coupler. This 
variation can be estimated by assuming the coupler in Figure 2.3 behaves as an acoustic 
waveguide where sound travels along the x direction as a plane wave and following the 
same procedure Zalesak uses to derive the wavelength dependent correction factor [9]. The 
acoustic pressure variation within the coupler is estimated as described earlier (Section 2.3), 
so Equations 2.10 and 2.11 still apply. Assuming U; = 0, the coefficient B was already 
determined to be 0. Now to solve for A, we again put 


Uz = —Su(L,t). (3.20) 


Applying this boundary condition to Equation 2.11, neglecting the time dependence term, 
and solving for A, we get 
—jwpUr pe 1 


a => j 21 
kSsin(KL) ~ S sin(kL)’ Cen 
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and Equation 2.10 becomes 


_ -jwpuy SiGe pe cos(kx) 
P= KS sin(kKL) ~ IS Sin(KL) 





(3.22) 


The expression for the pressure in the coupler at very low frequency (kL << 1), where the 
pressure does not vary over the length of the coupler can be derived from Equation 3.22 
using small angle approximations for the trigonometric functions, 

_ -jpcU, _ pe? 


Z ae 3.23 
PIE SkE jwV ee 





where the volume of the coupler, V, is equal to SL. Note that V/ pc? is the acoustic 

compliance of the coupler. Since the measured pressures are affected by this small variation 

across the coupler, they must be converted to low-frequency pressures to preserve the 

constant pressure assumption required in Equation 2.18 by multiplying by the ratio pj ¢/p. 

Dividing Equation 3.23 by Equation 3.22, we get this ratio, the correction factor F: 
Pif sin(kL) 


F=-1= 


p  kLcos(kx)’ G24) 


Zalesak [9] adds one additional refinement to the correction. In the actual coupler, the 
transducers do not form the walls at each end of the coupler, but are offset a small distance 
inward. He accounts for this difference by assuming that the coupler reciprocal transducers 
are infinitesimally thin and span the coupler cross section at a distance d from each wall of 
the coupler. Figure 3.5 shows these assumptions with the projecting transducers located at 
x = L—d and the receiving transducer at x = d. Applying the same low-frequency pressure 
assumption and ratio gives his wavelength dependence correction factor, Fy, given by 


sin(kL) 


~ KL cos(kd) cos(kx) ver 


Fy 
The additional cos(kd) term in the denominator improves the accuracy of F>. F and F are 
within 10% of each other over the H48 hydrophone frequency range 1-2000 Hz. Zalesak 
does not add any more features to the wavelength dependent correction analysis, although 
Figure 3.5 is still missing some features of the coupler. The reciprocal transducers located a 
distance d from each wall are actually 2.54 cm diameter spheres and the coupler cross section 
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is not constant — there is a small perpendicular leg in the middle, forming a T-shape where 
the H48 hydrophone is installed. The transducers are also assumed to be stationary. Zalesak 
argues that the correction is already small — for the H48 coupler, its maximum value is 
0.06 dB at 2000 Hz, falling to 0.01 dB at 700 Hz — and that further refinements would 
only make relatively small adjustments to the correction and are therefore negligible [9]. 
But additional refinements could be made using a lumped parameter, or more accurately, a 


finite element model. 





+X 


x=0 


Figure 3.5. Illustration of Zalesak’s coupler configuration assumptions for 
the wavelength dependent correction factor 


The ratio calculated to correct the measured pressures also applies to the measured voltage 
magnitudes. Two voltage measurements must be corrected: the hydrophone voltage and 
the receiving transducer voltage. The hydrophone voltage correction is given by Equation 
3.25 using x = L/2: 


sin(kL) _ 2sin (4) 


Dee a EEN, (3.26) 
kL cos(kd) cos (4) kL cos kd 


The receiving transducer voltage correction is also given by Equation 3.25 using x = d: 


sin(kL) 


Foy = ————_.. 
a kL cos?(kd) 


(3.27) 
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3.4 H48 Hydrophone Low-Frequency Equivalent-Circuit 
Model 


The behavior of the H48 hydrophone at low frequency, far below its first resonance, can 
be approximated as an RC high-pass filter circuit [13], shown in Figure 3.6. V is the open 
circuit voltage generated at the hydrophone, C the hydrophone capacitance, R the input 
resistance to the preamplifier, G the gain of the preamplifier, and V,,,;, the voltage output of 


the preamplifier. The response of the RC circuit is given by 


— Vout _ VxG 








M Se 3.28 
V 1+ 1/jwRC a 
and the phase angle by 
af 1 -1 ( fo 
1 1 
= ——]}= —|]. a.29 
é = tan (=z) tan (2] (3.29) 
The cutoff frequency, fo, is given by 
fo = (3.30) 
9 IARC’ 





Figure 3.6. Simplified hydrophone circuit model. Adapted from [13]. 


At the cutoff frequency, the sensitivity magnitude is -3 dB from its flat, high-frequency value 
and the phase is 45°. The low-frequency roll-off of the circuit is shown in Figure 3.7. Note 
that the magnitude response flattens out about one decade above the cutoff frequency, but 
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the phase does not flatten out until about two decades above the cutoff frequency. The low- 
frequency cutoff or 3 dB-down frequency of the H48 reference hydrophone was designed 
to be 0.035 Hz, but, after construction, was measured to be 0.08 Hz [14]. This (0.08 Hz) 
cutoff frequency implies that the phase response would begin to fall off below 8 Hz. At 1 
Hz, the phase should be about 18°. 
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Figure 3.7. Magnitude and phase roll-off at low frequency. Adapted from [15]. 


The simple RC high-pass filter model was fit to the complex sensitivity data using a 
nonlinear least-squares curve fitting procedure. For a given data set of input frequencies 
and complex sensitivity data, two fitting parameters were calculated from Equation 3.28, 
the voltage/gain product, V x G, and the RC time constant, RC. An unweighted, complex 
Levenberg-Marquardt algorithm was used, which, as implemented in Matlab’s /sqnonlin 
function [16], allows the least squares coefficients to be fit in the complex plane without 
having to separate their real and imaginary parts in the objective function. To allow for 
resistive and capacitive losses, the RC parameter was allowed to be complex, but the gain 
term was restricted to be a real number, as more detailed treatment of the preamplifier is left 
to future work. The fitted parameters were used to estimate the complex sensitivity, from 
which the the magnitude and phase were estimated by taking the respective magnitude and 
angle of the real and imaginary parts of the sensitivity estimate. 


Zi 


3.5 Procedure 

A complex-valued H48 coupler reciprocity calibration was performed across frequency, 
temperature, and pressure in accordance with [17]. Pressure was varied from 50 to 2000 
psi. Temperature was varied from 0°C to 40°C. For this report, at each temperature and 
pressure, a logarithmic frequency sweep was performed from 1-2000 Hz. At each frequency, 
20 cycles of the time series signals were recorded at a sample rate of five times the signal 
frequency. Each frequency sweep was repeated ten times. 
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CHAPTER 4: 
Uncertainty Estimation 





The new complex coupler reciprocity calibration procedure must account for the uncer- 
tainties treated by Zalesak [9] along with uncertainty in the phase response of the coupler 
components, the electronic equipment in the signal paths, the NI-6259 data acquisition 
hardware, and the effect of random error in calculating the complex sensitivity. The un- 
certainty is split into the two types described in the National Institute of Standards and 
Technology (NIST) guidelines for evaluating uncertainty [18], A and B. NIST type A un- 
certainty is calculated by performing statistical calculations on measurements. NIST type 
B covers uncertainties calculated in any other way, such as factoring in the accuracy of 


equipment used to make measurements or parameters used to calculate results. 


4.1 Estimation of NIST Type A Uncertainty 


Uncertainty in the complex sensitivity produces magnitude and phase uncertainty. Unlike 
measurements that fall into the set of real numbers, magnitudes must be greater than zero 
and phases fall in the range of - to 2, which is problematic for the arithmetic involved in 
uncertainty calculations. Analyzing the uncertainty in the complex plane, where the real 
and imaginary axes extend to positive and negative infinity avoids these problems [19]. 
Only after the uncertainties of the real and imaginary parts of the complex sensitivity are 
calculated are they converted to magnitude and phase uncertainties. 


The random error in the complex sensitivity is estimated using a bivariate normal distribution 
that accounts for the uncertainty in each of the real and imaginary parts as well as the 
correlation between the two, represented by a covariance matrix. Using singular value 
decomposition, the eigenvalues and eigenvectors of the covariance matrix can be used to 
define an elliptical region of uncertainty in the complex plane [19]. An example from an 
H48 data set is shown in Figure 4.1. The plot shows the real and imaginary parts of ten 
measurements of the complex sensitivity recorded at a single frequency. This data group 
and confidence region represent one of the data groups and confidence regions found in 
Figure 5.1, a data set recorded across the frequency range 5 Hz - 2 kHz. The lengths of the 
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major and minor semiaxes are related to the eigenvalues of the covariance matrix, and the 
rotation angle is determined by the correlation of the real and imaginary parts. The length 
of the major semiaxis, a, is given by the largest eigenvalue and the minor semiaxis, b, by 


the smallest, as shown in Equations 4.1 and 4.2: 


a= WVSAmax 5 (4.1) 
b=WVsAmin 5 (4.2) 


where s is a coverage factor given by 


we ny 212), (4.3) 
determined by the F>,-2 inverse cumulative distribution function with two degrees of 
freedom and n observations at a level of confidence a [20]. The F,-2 is used instead of 
the standard normal distribution used for scalar measurements, resulting in coverage factors 
different than those used for scalar measurements, (i.e., k = 1.96 for 95% confidence). The 
ellipse in Figure 4.1 is scaled using Equation 4.3 with a 95% confidence interval resulting 
in a coverage factor of 2.45 for 1000 observations. The major axis of the ellipse is rotated 
from the real axis by the angle between the largest eigenvector of the covariance matrix and 
the real axis [21]. 


The uncertainty in the magnitude and phase can be calculated from the magnitude and 
phase boundary points on the covariance ellipse shown in Figure 4.1. The magnitude 
boundary points are selected to give the longest and shortest phasors (from the origin to 
the covariance ellipse boundary) and the phase boundary points, to give the smallest and 
largest possible angles from a phasor on the covariance ellipse boundary to the real axis. 
These minimums and maximums were calculated using fminbnd in Matlab and are used to 
define the uncertainty at a particular frequency. The magnitude and phase errors across the 


frequency range are shown in Figure 4.2. 
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Figure 4.2. Magnitude and phase random errors across frequency 
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4.2 Estimation of NIST Type B Uncertainty 

Random error from the complex voltage measurements contributes to uncertainty in the 
magnitude and phase of the sensitivity, but other factors will contribute additional uncer- 
tainty. The total phase uncertainty is that calculated from the covariance matrix depicted 
in Figure 4.2, combined with the uncertainty in the phase measurements taken with the 
HP35565 signal analyzer, +0.5°, used to correct for the phase shift in the differential 


Obtor = 5b? nq + 0-52. (4.4) 


The final magnitude uncertainty is a combination of uncertainties calculated from the 


amplifier, 


covariance matrix and additional uncertainties from measured quantities in Equation 2.18, 


as shown in Equation 4.5 (complex values assumed here): 











2 2 
6|M, 6|M, 6|M, 
|My | ( | a) (4 a) (4.5) 
cov add 


[Mal tora = V\ Ma] [Mu 


The additional uncertainty term can be broken into electrical and acoustic compliance 


groupings, with an uncertainty term assigned for each parameter in Equation 2.18. 


In the electrical group, the uncertainties in voltage measurements from the hydrophone, 
deny and deny arise from the accuracy limitation of the NI-6259 used to make the voltage 
measurements. The uncertainties in voltage measurements from the receiving reciprocal 
transducers, dexy and deyy, arise from the accuracy limitation of the NI-6259 along with 
the gain accuracy of the Ithaco 1201 preamplifier used to amplify the signal in front of the NI- 
6259. The uncertainties in voltage measurements across the standard capacitor, decy and 
décy, also arise from the accuracy limitation of the NI-6259 with an additional contribution 
from the differential amplifier used amplify those signals. Finally, the uncertainty in the 
capacitance of the standard capacitor as manufactured, 6C, must be accounted for. These 


terms are explained in more detail in Section 4.2.1. 


In the acoustic compliance group, there is uncertainty in the volume of the coupler cavity, 
OV, associated with the method used to measure that volume. The uncertainty in the speed 
of sound in castor oil, 6c, arises from uncertainties in the estimate of that value [22] and in 


the uncertainty in the temperature and pressure inside the coupler cavity. The uncertainty in 
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Figure 4.3. Components of fractional uncertainty in the sensitivity magnitude 
from Equation 4.6, taken from a data set at 1000 psi and 3.2°C 


the density of castor oil, 6, also arises from the uncertainty in the temperature and pressure 


inside the coupler cavity. These terms are explained in more detail in Section 4.2.2. 


The size of the fractional uncertainty magnitudes of each of these additional uncertainties 
as components of the sensitivity magnitude uncertainty are shown in Figure 4.3, based on 
the relationships of the quantities in Equation 2.18. The fractional uncertainty in the sound 
speed dominates and can be as high as 3% at the top of the temperature range (note that the 
effect of the fractional uncertainty in the sound speed is doubled). The volume fractional 
uncertainty is less than half the amount of the sound speed. The fractional uncertainties in the 
voltages measured across the standard capacitor are a close third, despite their contribution 
being reduced by half. Those voltage uncertainties are dominated by the 1.6% measurement 
error in the differential amplifier. Similarly the next largest fractional uncertainty, that of the 
receiving reciprocal transducer measurement, is dominated by the 1% error in the gain of 
the Ithaco 1202 preamplifier. The fractional uncertainty in the capacitance of the standard 
capacitor used in the current measurement is given by the manufacturer to be 0.1%. The 
smallest contribution comes from the fractional uncertainty of the hydrophone voltages 
measured by the NI-6258. 
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These additional fractional uncertainties can be combined in quadrature for a total additional 


fractional uncertainty of 1.4%, as shown in Equation 4.6: 
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4.2.1 Uncertainty in the measured voltage magnitudes 





The magnitude of the voltage terms eyx, eny, exy, cx, eyx, and écy all contain uncertainty 
due to the absolute accuracy of the NI-6259, which varies with temperature and can be 
calculated from the NI-6259 specifications [10] as shown in Equation 4.7: 


Accuracyny.6259 = Vmeas(gainErr) + + (offsetErr) + noiseErr. (4.7) 


For the coupler measurement setup, this uncertainty is split roughly equally into an offset 
error, OffsetErr, and gain error, gainErr, with gain error being slightly smaller at room 
temperature. The constant 1/2 is selected from a chart of values in [10] and changes based 
on the selected input range in the NI-6259. For the number of readings averaged for each 
voltage measurement, the noise error term is negligible. Offset and gain errors are shown 
in Equations 4.8 and 4.9. 

ppm 


°C 


ppm 


gainErr = 80ppm + (13 x TempDrift) + (1 OG x CalTempDitt] (4.8) 


offsetErr = 20ppm + (27 x TempDrift) + 60ppm (4.9) 


TempDrift is the change in temperature since the last internal calibration. The NI-6259 
can perform an internal calibration on command that resets this temperature difference, 
which can be performed at the start of each day’s measurements. Therefore TempDrift 
amounts to the change in temperature in the coupler lab over the course of the day’s 
measurements. CalTempDiff is the difference in temperature between the coupler lab and 
the lab at National Instruments where the NI-6259 was originally calibrated. Note that the 
constants are selected from a chart of values in [10] and changed based on the selected input 
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range in the NI-6259. The offset portion of the error depends on the selected input range (a 
constant) and temperature of the NI-6259. Over the short time period it takes to collect a 
set of measurements at one frequency, this error will tend to shift all the measurements in 
the same direction. The least-squares fit on that frequency data set will separate that shift 
into a DC offset term, which is discarded from sensitivity calculations. Therefore the offset 
portion of the NI-6259 uncertainty is ignored. The remaining gain portion will distort only 
the magnitude of the signal and will be factored in to the voltage uncertainties in Equation 
4.6. 


The voltage magnitudes representing the current to transducers X and Y (|ecx| and |ecy|) 
have an additional uncertainty due to the differential amplifier used to make those measure- 
ments. These voltage measurements are made across a standard | uF capacitor through a 
differential amplifier, as shown in the simplified circuit in Figure 4.4. 


Vs 


Vout 





Figure 4.4. Simplified circuit of differential amplifier, standard capacitor, 
and projector 


Given a source voltage, V;, and a capacitance of 4.32 nF for the projector (Cp;o;) and 1uF 
for the standard capacitor (Cs;q), the voltage, vo,;, measured by the differential amplifier is 
given in [23] as 

Vout = Ad(U1 — 02) + AcmUem > (4.10) 


where Aq is the differential gain, A,,, is the common-mode gain, and v,,, is the voltage 
common to both inputs of the amplifier, which in this case, is just v2. From the capacitances, 


the voltage drop across the standard capacitor, vj — v2, is 0.0043V;, and the voltage drop 
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across the projector, v2, is 0.9957V;. The error term is the common-mode contribution to 


the voltage output, which can be expressed as the fractional error 





Acm Vem 
Verr = — —— . 4.11) 
Aa(v1 — v2) 
The common-mode rejection ratio of the differential amplifier is defined in [23] as 
Ad 
CMRR = 20 log Acnl (4.12) 


and is given by the manufacturer to be -83 dB. Solving Equation 4.12 for the ratio Ag/|Agm| 
and inserting into Equation 4.11, we get a voltage error of 1.6% [9]. This uncertainty is 
added in quadrature with the uncertainty from the NI-6259 for these two voltage terms. 


4.2.2 Uncertainty in the acoustic compliance terms 
The volume, density, and sound speed terms enter the sensitivity calculation by way of the 
acoustic compliance term in the reciprocity parameter, introduced in Equation 2.17 and 


reproduced in Equation 4.13: 
V 
peo’ 


where V is the volume of the castor oil filling the coupler cavity, p is the density of the 


Ca = (4.13) 


castor oil, and c is the speed of sound castor oil. 


The fractional uncertainty in the volume is 1% [9] based on the accuracy of the original 


volume measurement. 


The sound speed in castor oil is estimated from the empirical relation in Equation 4.14 


given by Timme [22], with an absolute accuracy of +1.1%: 
c(T, p) = 1570(ap + ayT + axl +a,Tp+aap + asp’). (4.14) 


Inputs to this relation are temperature, T, in degrees Celsius, pressure, p, in MPa, and the 
constants dg through as given in Table 4.1. The ambient temperature and pressure are only 


held constant within +1°and +10 psi, respectively, contributing uncertainty to the calculated 
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sound speed as shown in Equation 4.15: 
5c(T, p) = 1570 ((a; + 3a3T? + ap) ST + (ayT + ag + 2asp)dp) . (4.15) 


These errors do not add in quadrature as they are not independent. Changes in temperature 
of the coupler are directly proportional to changes in pressure within a data set. The absolute 
accuracy of the sound speed equation is added in quadrature to the fractional uncertainty 
due to temperature and pressure variation for the total fractional uncertainty in the sound 


speed shown in Equation 4.16: 


Cc 


2 
6 (AS) + (1.1%)2. (4.16) 


Table 4.1. Constants for the empirical relationship of sound speed in castor 
oil with temperature and pressure 











Constant Value Units 
ao 1 none 
4 1 
ay ~2.15x 10 eee 
°C 
4.0 x 107° : 
a2 OX oC? 
a 2.5 x 10~° : 
? °C x MPa 
1 
2.22 x 1077 
oe MPa 
1 
-~3.0 x 10-6 a 
2 MPa 





a1 


The density of castor oil is estimated by the relation given in Stallard [24] as 


1 
— = by + boT + b3T? + baT? + bsp + bopT + b7pT? + bepT?+ 
p 
bop” + biop°T + bup’T? + dip? + bisp°T, (4.17) 
where T is the temperature in degrees Celsius and p is the pressure in MPa. The coefficients 


are given in Table 4.2. As with the sound speed estimate above, uncertainty is added due to 


the variation of temperature and pressure within the test setup, given in Equation 4.18: 


dp = 64 = (by + 2b3T + 3b4T? + bop + 2b7pT + 3bgpT? + biop” + 2bip°T + bi3p*) ST + 
(bs + beT + b7T? + bgT? + 2bop + 2byopT + 2b\;pT* + 3bi2p” + 3b13p'T) dp. (4.18) 


These errors do not add in quadrature, as changes in temperature of the coupler are directly 


proportional to changes in pressure. 


The uncertainty in the compliance term can vary with pressure and temperature. The varia- 
tion in uncertainty across the pressure range is negligible. The variation across temperature 


is shown in Figure 4.5. 
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Table 4.2. Constants for the empirical relationship of castor oil density with 


temperature and pressure 








Constant Value Units 

by 1.027 cneg7! 

bo 7.038 x 1074 cnbg (°C)! 

b3 9.659 x 1077 cm3g7! (°C)? 

ba 3.045 x 107? cme g (°C)? 

bs —4.909 x 1074 cng !M Pa"! 
be —~2.633 x 10° = em?g™!MPa!(°C)"! 
by —4.042x 107? = em3g™!MPa7!(°C)~? 
bg ~8.772x 107!  em3g7!MPa7!(°C)-3 
bo 1.471 x 107° cm? g!MPa~ 
bio 9.181x 107°? cm3g7!MPa?(°C)! 
bi 3.591107! = em3g™!'MPa!(°C)~? 
by —3.634 x 107° cn?g-!MPa~> 

by3 —1.657x 107! em3g™!MPa3(°C)"! 
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Figure 4.5. Fractional uncertainty in the components of the acoustic 
compliance in the calculation of sensitivity 
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CHAPTER 5: 
Results and Discussion 





The results of two complex hydrophone calibrations with systematic magnitude and phase 
corrections, NIST type A and B uncertainties, and fitted simple, high-pass equivalent- 
circuit models are presented. The complex sensitivity is plotted in the complex plane, and 


the magnitude and phase is plotted with respected to frequency. 


The complex free-field voltage sensitivity of H48 serial number two with the associated 
covariance ellipses and the fitted simple, high-pass equivalent-circuit model is shown in 
Figure 5.1. Each ellipse represents a confidence region around the average sensitivity 
measurement representing the NIST type A error at a particular frequency between 5 Hz 
and 2000 Hz (spaced logarithmically), calculated as shown in Section 4.1. Each average 
measurement is the mean value of ten separate fittings of the least squares parameters to a 
sampled waveform. Each least squares fit is performed on data sampled at five times the 
signal frequency for 20 signal periods, or 100 data points. The frequencies begin at 5 Hz 
at the smallest values of the real and imaginary parts of the sensitivity in Figure 5.1. As 
frequency increases, both the real and imaginary parts of the sensitivity increase. As the 
imaginary part of the sensitivity approaches zero above 100 Hz, the real part decreases until 
the frequency reaches 2000 Hz. The model is a good fit for the data at low frequency, but 
fits poorly above 100 Hz, above those of primary interest in this report. The relationship 
between the real and imaginary parts of the complex, open-circuit voltage sensitivity are 
shown in Figures 5.2 and 5.3, respectively, where the frequency above which the model fit 
deteriorates is more apparent. More physically intuitive conclusions can be drawn from 


plotting the magnitude and phase of the data and model. 


The magnitude portion of the free-field voltage sensitivity (FFVS) is calculated by taking the 
magnitude of average sensitivity values from Figure 5.1 and converting to dB re 1V/uPa. 
The FFVS magnitude of H48 serial number 4 is shown in Figure 5.4. The magnitude of the 
model is calculated in the same way from the estimated complex FFVS and converted to 
dB re 1V/uPa. These data can also be verified by comparison with the FFVS magnitude 
calculated from RMS voltage measurements made with the multimeters in a standard 


4] 


imaginary 


Confidence Region 
* Average Value 
——r= Mode! 
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real ~ 104 


Figure 5.1. Complex open-circuit voltage sensitivity measurements and un- 
certainty fitted with a simple circuit model for H48, serial number two, from 
a data set taken at 500 psi and 39.6°C 


coupler calibration [17]. The FFVS magnitudes measured from each method agree within 
measurement uncertainty above 15 Hz. Below 15 Hz, the signal is below the minimum 
frequency the HP 3478A multimeters are certified to measure, which may explain the 
disagreement in the data and the upward trend in the FFVS calculated from multimeter data 
at low frequency. This upward trend is not supported by the model nor by the H48 design 
document that specifies the typical FFVS for an H48 hydrophone [14]. It is not clear from 
the sensitivity magnitude data where the cutoff frequency is for the hydrophone. At the high 
end of the frequency range, both sets of measurements show a pronounced hump where the 
model is flat. This poor fit corresponds to the hooked section of the model in Figure 5.1. 


The phase response and model are calculated by taking the angle of the complex values and 
their modeled estimates in Figure 5.1. The phase response of H48 serial number four is 


compared to the simple, high-pass equivalent-circuit model in Figure 5.5. At low and high 
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Figure 5.2. The real part of the complex open-circuit voltage sensitivity 
measurements plotted versus frequency with NIST type A uncertainty and 
the simple, high-pass equivalent-circuit model for H48, serial number 2, from 
a data set taken at 500 psi and 39.6°C 


frequencies, the model is in good agreement with the data. In the middle of the frequency 
range, the phase data does not flatten out like the model. The fit in this range can be 
improved by letting the gain coefficient of the simple circuit model be complex, introducing 
a stray capacitance-like term to the input of preamplifier. This analysis is excluded, as a 


more detailed model of the preamplifier is left to future work. 


At 1 Hz, the measured and modeled phase is 20°. It should be noted for low frequency 
calibrations that the constant phase assumption breaks down between 10 Hz and 20 Hz 
and gets increasingly worse as frequency continues to decrease. Since the phase begins its 
roll-off at a higher frequency, it is easier to calculate the cutoff frequency with these data 
than with the magnitude data alone. Extrapolating with the model to a phase of 45°, we get 
an estimated cutoff frequency of 0.3 Hz, higher than the 0.08 Hz measured when the H48 
hydrophones were built [14]. Assuming even 10° of error in the model extrapolated down 
to 0.3 Hz gives an uncertainty in the estimated cutoff frequency of 0.1 Hz, as seen in Figure 
3.7. 
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Figure 5.3. The imaginary part of the complex open-circuit voltage sensitivity 
measurements plotted versus frequency with NIST type A uncertainty and 
the simple, high-pass equivalent-circuit model for H48, serial number two, 
from a data set taken at 500 psi and 39.6°C 
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Figure 5.4. Open-circuit voltage sensitivity magnitude of H48, serial number 
four, comparison between times series measurements, RMS measurements, 
and a simple circuit model for a data set at 1000 psi and 3.2°C 
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Figure 5.5. Phase response of H48, serial number four, comparison between 


times series measurements and a simple circuit model for a data set at 1000 
psi and 3.2°C 
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CHAPTER 6: 


Conclusion 





A determination of the phase as well as the magnitude of a hydrophone’s free-field voltage 
sensitivity was incorporated into the USRD reciprocity coupler calibration procedure by 
using the hydrophone sensitivity equation for a reciprocity calibration, as shown in Equation 
2.6, with complex voltage inputs and a complex reciprocity parameter, given in Equation 
2.17 . The complex input signals are calculated from time series measurements using a 


least squares estimator. 


The physical setup remains the same as the previous, magnitude only, calibration setup, 
with data acquisition equipment added to record time series data in place of the RMS 
multimeters. The electronics add their own phase errors to the system, which must be 
separately measured and removed to isolate the phase response of the hydrophone under 
test. 


Data sets were recorded back-to-back on H48 hydrophones using the standard RMS mul- 
timeters and new data acquisition equipment to record time series measurements. The 
magnitude calculated from these data sets agreed within measurement uncertainty. Both 
magnitude and newly incorporated phase data were fitted to a simple, RC high-pass filter 
equivalent circuit model of the hydrophone behavior at low frequencies. The magnitude 
data fit the model within measurement uncertainty across the applicable frequency range, 
although the hump in the higher frequencies was not predicted. The phase data fit well in the 
lower frequency range, of greatest interest in this report, but did not fit within measurement 


uncertainty in the middle of the frequency range. 


Both magnitude and phase response began to roll off as frequencies approached | Hz from 
above, but the trend is easier to see in the phase data, which begins to roll off a decade higher 
than the magnitude data for a given cutoff frequency. The phase data were extrapolated 
with the equivalent-circuit model to estimate the cutoff frequency of the hydrophone — 
0.3 Hz for H48 SN4 — which is greater than measured previously when the H48 was first 
constructed [14]. The higher cut-off frequency could lead to inaccuracies when measuring 


phase with an H48 hydrophone or associated secondary references under the assumption 
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that phase is constant down to 1 Hz. 


Errors and uncertainties in the calibration are similar to that calculated previously [9] for the 
magnitude-only system, with new errors for the new data acquisition equipment. Additional 
errors and uncertainties had to be accounted for in the phase measurements, the majority of 
which come from the electronic equipment. Uncertainty from random error was calculated 


from the covariance matrix of the complex sensitivity data. 


Additional acoustic sources of error come into play higher in the H48 frequency range, and 
should be modeled in future studies of this coupler. The irregularities in the shape of the 
coupler, the size of the reciprocal transducers compared to the coupler and their position 
within, and viscous effects can be estimated with a lumped parameter model, but would be 


best addressed through finite element analysis. 
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CHAP TER 7: 
Future Work 





Some features of the coupler do not conform to the simplifying assumptions of the coupler 
diagram in Figure 2.3 and lead to errors not treated in this report. These errors arise from 
the fact that the reciprocal transducers are 2.54-cm diameter spheres mounted inside the 
coupler walls, and the coupler itself has a shallow T-shaped branch at its center, where 
the hydrophone mounts. Future work should address these differences through a lumped 
parameter, or preferably, a finite-element model. Figure 7.1 shows a scale diagram of a 
cross section through the center of the coupler cavity. X and Y are the spherical elements 


of the reciprocal transducers and H is the hydrophone element. 


The spherical transducers mounted inside the coupler walls are not necessarily stationary 
and do not occupy the full cross section of the coupler, leaving a small cavity between them 
and each end of the coupler. Incomplete cross sectional coverage leads to a radiation mass 
term in the acoustical transfer impedance. Space on either side of the reciprocal transducers 
where the flow converges and diverges contributes inertia and a Helmholtz resonance. 
Another divergence in the flow occurs at the T-shaped branch. These contributions can 
be captured by modifying the coupler compliance to be a frequency dependent, effective 
compliance. These effects only contribute to the acoustic radiation reactance of the coupler 


and therefore only affect the magnitude of the coupler calibrations, not the phase. However, 


Figure 7.1. Coupler scaled diagram with reciprocal transducers X and Y, and 
hydrophone H 
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viscous effects not presently accounted for in modeling the acoustic field within the coupler 


could introduce phase errors. 


Shortcomings in the simplified circuit model probably lead to a poor fit in the phase data 
in the middle of the frequency range. Adding more detail to the model of the preamplifier 


could help explain these discrepancies. 


Also worth noting, in cases where the NIST type B error is much larger than the NIST type 
A error, fitting a model using an unweighted, nonlinear least-squares approach, as done 
in this report, will give similar results as a weighted, nonlinear least-squares fit. In cases 
where the NIST type A error is more significant, using a weighted fitting algorithm should 


improve the model results. 
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